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Abstract. - Soft interfaces can mediate interactions between particles bound to them. The 
force transmitted through the surface geometry on a particle may be expressed as a closed line 
integral of the surface stress tensor around that particle. This contour may be deformed to 
exploit the symmetries present; for two identical particles, one obtains an exact expression for 
the force between them in terms of the local surface geometry of their mid-plane; in the case of 
a fluid membrane the sign of the interaction is often evident. The approach, by construction, 
is adapted directly to the surface and is independent of its parameterization. Furthermore, it 
is applicable for arbitrarily large deformations; in particular, it remains valid beyond the linear 
small-gradient regime. 



Introduction. - Physical interactions between spatially separated particles are mediated 
by fields: matter interacts by curving spacetime itself; if it is charged it will interact through 
the electromagnetic field. Equally well, interactions of an indirect nature may be mediated 
by effective fields: for example, colloidal particles in suspension can interact by disturbing 
local ion densities (if they are charged) [1], by distorting the order of an embedding liquid 
crystal [2] , or by locally phase separating a binary mixture [3] . 

An important class of such interactions, purely geometrical in nature, are those between 
particles localized at an interface which are induced through the local deformation in its shape. 
The simplest example is provided by capillary interactions, originating in surface tension with 
an energy proportional to the excess area of the interface. These interactions play a key role 
in many technological processes, among them ore fiotation and foam stabilization [4]. To 
describe the compression and bending of surfactant monolayers and lipid membranes [5-8], 
surface energies involving higher order derivatives are required. For example, the interactions 
between membrane bound proteins, which play an important role in cell biology [9], are 
described by an energy quadratic in the curvature of the membrane surface. 

A major obstacle to providing a theoretical description of these surface mediated inter- 
actions is that the corresponding field equations are nonlinear. It is easy to see why: while 
the underlying surface free energy density may be simple — some quadratic invariant of gen- 
eralized strains and thus "harmonic" — the description of the curved manifold it inhabits is 
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intrinsically non-linear. Even in the case of an essentially flat surface, where its description in 
terms of a height function ( "Monge gauge" ) keeping only the lowest order in gradients leads 
to approximate linear differential equations (and thus to Green functions), exact analytical 
solutions are difficult to obtain. This is because a straightforward superposition ansatz [10] 
will typically violate the correct boundary conditions. It remains entirely unclear how any 
linear result generalizes to the full nonlinear situation. 

Despite this obstacle, substantial progress can still be made by recognizing that the deter- 
mination of the interaction between two objects mediated by the surface involves solving two 
distinct problems: first, the particular way in which objects bind to the interface, together 
with its elastic specifications, determines its equilibrium {i. e. energy minimizing) shape — find 
that shape; second, equilibrium of the resulting complex generally requires additional external 
forces which constrain the bound objects at their positions — find these forces. Usually one 
thinks of the second problem in a way which depends heavily on the successful implementation 
of the first: since the energy of the complex depends on the relative positions of the bound 
objects, appropriate derivatives of the energy with respect to their coordinates yield the forces 
one is looking for. Rarely, however, is it possible to write down this "potential" energy; as a 
result, this conceptually straightforward route to the forces is infeasible in practice. In other 
words, a solution to part 1 is generally impossible, leaving no further handle on part 2. This 
reasoning distracts from the fact that forces transmitted by an interface must be encoded 
directly in its shape; if the Hamiltonian describing the surface is geometrical, so also are the 
stresses underpinning this shape. 

In this letter we reformulate the full nonlinear interaction problem in terms of surface 
geometry; the conserved covariant surface stress tensor [11] associated with the geometry will 
play a central role. We show how to express the surface mediated force on a particle as a 
closed line integral of the stress tensor. By suitably deforming the contour to exploit the 
symmetries of the configuration, a remarkably simple and transparent expression for the force 
is obtained. In particular, the sign of the force may turn out to be evident using only very 
qualitative features of the geometry. We should mention that a stress tensor approach to 
membrane mediated interactions was suggested by Kralchevsky et al. [6]. However, the full 
potential of a completely geometric description has not previously been exploited. 

Differential geometry, the energy, and the stress tensor. - In this section we sketch the 
essential geometric background. While this will define our notation, it is most likely too 
concise to serve as a stand-alone introduction. For the geometry we therefore refer the reader 
to Refs. [12]; the stress tensor is introduced at greater length in Ref. [11, 13]. 

Consider a surface S embedded in three-dimensional Euclidean space M'', which is de- 
scribed locally by its position X(^^,^^) G R^, where the are a suitable set of (curvilinear) 
local coordinates on the surface. Given the two tangent vectors = dX/d^°' = daX and 
the (unit) normal vector n = ei x e2/|ei x 62!, the surface geometry is described completely 
in terms of the induced metric gab = Ba ■ and the extrinsic curvature Kab = —n ■ da^b [14]. 
We denote the metric-compatible covariant derivative by Va and the corresponding Laplacian 
as A = VoV". As usual, repeated indices - one up one down - imply a summation. The total 
curvature K is the trace of the extrinsic curvature, K = g°''Kab [15]. 

We now associate with the surface an energy which can be written as a surface integral 
over a scalar Hamiltonian density Jf constructed out of local geometric invariants. 




(1) 



where the infinitesimal surface element is dA = d^^^/g with g = det{gab)- Let us perform a 
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Table I - Euler-Lagrange derivative £'{Jff) and components of the stress tensor /" — f°'^eh~\- f^n for 
several simple scalar surface Hamiltonian densities . Projecting Vaf"' ~ Sn implies the identities 
Va/" - Kahf'' = and Va/"'' - Klf" =0. R is the intrinsic scalar curvature [15]. Notice that 
and K°'^Kab yield identical S and /" (a consequence of the Gauss-Bonnet-Theorem [12]). These 
results may be derived using techniques developed in Refs. [11,13]. 



variation X X + 6X of the embedding functions. The concomitant first variation of the 
functional ^ can be cast as a bulk part plus a pure divergence: 

5Hs ^ [ dA S{.J^) n-5X + [ dA VaQ" ■ (2) 

Here, S'n is the bulk Euler-Lagrange derivative of i which is evidently purely normal. The 
second integral (which is identical to the boundary integral of over 9E) originates in the 
tangential variations as well as the derivatives of normal variations: 

(5° = — /° • SX + possible terms containing derivatives of SX . (3) 

The object is the surface stress tensor. Its components in the local frame {ei,e2,n}, as 
well as the Euler Lagrange derivative S'{ J^f), are listed for a few simple Hamiltonian densities 
in Table n 

Suppose SX is simply a constant translation, which of course leaves the Hamiltonian 
invariant. We thus have 5H^ = 0, and with the help of Eqns. (03) Vq/" = S'{Jf)n. But 
a true equilibrium surface is stationary with respect to arbitrary variations. Thus the Euler- 
Lagrange ("shape") equation S' — Q also holds, and we get the conservation law 

= . (4) 

Its existence is simply a consequence of Noether's theorem: a continuous symmetry implies 
a conservation law on shell. Note that if the surface encloses a fixed volume V , a term 
—PV = —\P dA n ■ X involving the Lagrange multiplier P needs to be included in the 
functional, yielding the shape equation S' = P, and Eqn. Q is replaced by V af^ — Pn. The 
same holds if there exists a pressure drop P across the two sides of the interface. 

Forces via the stress tensor. - In elasticity theory the divergence of the stress tensor 
equals the external force per unit volume of a strained material [16]. Likewise, the divergence 
of the surface stress tensor defined above equals the force per unit area of the strained surface. 
For instance, the equation Va/° ~ Pn then states that a pressure P across the surface is a 
source of stress. Using Stokes' theorem, the total force acting "within" any patch Sq is 



Fso ^ f dA Vaf ^ i ds laf 



(5) 
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Fig. 1 - Illustration (view from the top) of how the force on one of a pair of objects bound to an 
interface can be calculated as a closed loop integral of the surface stress tensor. The contour can 
subsequently be deformed to conform to the symmetry of the situation. 

where I = I'^Ba is the outward pointing unit normal to the boundary curve 9Eo (which is 
by construction tangential to the surface), and s is the arc- length on dY,a. If So is a free 
equilibrium patch {i. e., no external stresses), Eqn. Q shows that F^g — 0. Generally, 
however, the patch will contain regions where external stresses act and i^s^ will be nonzero. 
Observe now that for P = Q the total force equals the line integral of the stress tensor along 
any curve enclosing these sources. This is because V af^ — permits us to deform the contour 
of integration — provided we do not cross any of these sources of stress in doing so. The case 
of P 7^ will be treated elsewhere. 

Consider now two identical particles bound to some asymptotically flat surface, as is 
schematically sketched in Fig. ^ Due to surface mediated interactions such a situation can 
only be stationary if external constraining forces fix the particle positions. These forces trans- 
mit stresses onto the surface, which are thus picked up by a line integral over the surface stress 
tensor around either particle. In the absence of a pressure difference the contour of the line 
integral can be deformed so as to take advantage of the available symmetry, see again Fig. ^ 
Once the contour is pulled open wide enough, the surface will be flat at branches 2, 3, and 
4, and the stress tensor will be very simple. In fact, the contributions from branch 2 and 4 
will then cancel each other. The only nontrivial contribution stems from branch 1, and its 
evaluation is greatly simplified by the symmetry. 

This outlines the basic strategy, which can evidently be tailored towards many other 
situations. We will now demonstrate its application to a few important standard cases. 

Force in fluid membranes. - The example of a surface we would like to focus on is an 
elastic symmetric fiuid membrane, described by the surface Hamiltonian [17] 



Here, k is the bending stiffness and a the surface tension. For the special case k — this 
reduces to the problem of a surface with surface tension only, and describes a soap film on 
large enough length scales or a water surface on length scales smaller than the capillary length. 
From Table CI we find that the associated stress tensor is given by 




(6) 
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Fig. 2 - Illustration of the geometry of a symmetric (solid line) and an antisymmetric (dotted line) 
two-particle attachment. The angle between surface normal at the symmetry point and vertical is tp. 



We now introduce (orthonormal) tangent vectors {t,l} along branch 1: t points along the 
integration line (t — f^ea), and (recall) I points normally outward. A short calculation then 
shows that the force stemming from branch 1 is given by 



Fi = - ds 



I - K{W_LK)n \ , (8) 



where we have defined the principal curvatures K± = l°-fKab and K\\ — t'^t''Kab, as well as 
the derivative along I, V± = I'^Va = d/dl. By symmetry there is no contribution to Fi along 
t. The minus sign out front stems from the fact that the membrane mediated force on the 
particle is opposite to the external force necessary to counterbalance it. 

To proceed further, we must look separately at the two different possible symmetries: 
either two particles adhere at the same side of the membrane (symmetric) or at different sides 
(antisymmetric), see Fig. In the symmetric case the curvatures K± and both have an 
extremum in the I direction along branch 1, so W±K = and the normal force component 
vanishes there. Furthermore, since the profile is horizontal in the middle, I = x there (where x 
is the unit vector pointing in the horizontal x-direction) . Finally, on branch 3 the stress tensor 
simplifies to 3 = —aSa, since far enough away the curvature becomes zero. Therefore, the 
total force ^1+^3 = i^sym^^ on the left particle is given by 



Fsym = (tAL 



\Kjds{Kl-Kl), (9) 



where AL > is the excess length of branch 1 compared to branch 3. We immediately see that 
the contribution due to tension is attractive. The curvature contribution, on the other hand, is 
the integral over the difference between the squared principal curvatures along the mid-curve, 
and as such has no evident sign. However, if we had adsorbed two parallel cylinders, which 
are sufficiently long such as to neglect end effects, the contribution if^ vanishes. Furthermore, 
the mid-curve becomes a line, and thus AL = 0. In this case we find for the force per unit 
length of the cylinder 

Fsym.cyl/L = -^^Kl < , (10) 

which is evidently always repulsive. Note that even though the tension a does not occur 
explicitly, it will enter the force indirectly through its influence on the value of K± . In Monge 
gauge, with height function h, K± = [h'/{l + h''^)]'\x=o = h"{0), since h'{0) = 0. Using this, 
Eqn. (|10|l is then quantitatively corroborated (in the linear regime) by the calculations in 
Ref. [8]. Details of this will be presented elsewhere. 

In the antisymmetric case branch 1 is a rotational symmetry axis of degree 2 and thus 
a line; hence both Kn and K± vanish. Since K± changes sign from positive to negative, 
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VxXj^ < 0. Observe that the profile on the midUne is always tilted in the direction indicated 
in Fig. 121 because otherwise it would have three nodal points on the asymptotic horizontal line 
and not just one, and the energy is expected to be higher. If we conceive of the constraining 
external forces as fixing the horizontal separation but not the vertical position of the particles, 
the latter is equilibrated and the vertical force component vanishes. The force which remains 
on the left particle is thus again horizontal, .Fantisym = ^antisymaJ, and given by 

^^antisym ^ ^'^ Cr ( COS (^(s) - l) ~ K siu (^(s) (if ^ + ) . (11) 

This time the tension contribution is repulsive. In the bending contribution the V±K± is 
attractive, but unfortunately not much can be said about the sign of Vj_i4r||, hence the overall 
sign is not obvious. However, in the case of two cylinders bound on opposite sides, Eqn. Q 
implies that |/"| is constant on each of the three membrane segments, and the force expression 
(per unit length) simplifies to [18] 



Fantisym,cyl/i = I /midpoint I " ^ = - <J > 0, (12) 

which is manifestly positive, implying particle attraction. If we expand the square root to 
first order in the bending part, we arrive at the result for linear theory i^antisym.cyi,^ <c i/-^ = 
^k{W±K±)'^ (with A = y^K/a). This latter expression is again quantitatively confirmed by 
using profile and force as they are calculated in Ref. [8]. 

Discussion. - In the previous sections we have outlined a very general method to obtain 
exact results for surface mediated interactions which sidesteps the need to solve the field equa- 
tions explicitly. This approach is fully covariant: one is free to choose the parameterization 
which is most appropriate; it is also valid for large deformations and not limited to a linear 
approximation. To our knowledge, there are no analogues based on energy minimization of 
our results outside this regime. It is true that our formulas contain unknown quantities re- 
lated to the shape of the surface. It was clear from the beginning, however, that we could 
not expect to solve the problem completely without determining this shape; what has been 
established is the connection between the geometry of the surface and the forces transmitted 
by it, something which does not come across when energy is differentiated - even in the linear 
regime. With sufficient ingenuity, it may be possible to extract the necessary information on 
the shape from the Euler-Lagrange equation (which so far we have not used at all!) without 
needing to solve it explicitly. Even without this input the functional form alone may identify 
the sign of the interaction [19]. Most importantly, this framework provides a new window 
onto an old and important problem. It can be combined with any approach, be it analytical 
or numerical, which determines the surface shape. At the very least, we have shown that it 
provides valuable non-trivial consistency conditions for analytical calculations. When such 
calculations are ruled out, the determination of the force between particles by numerically in- 
tegrating the stress around one of them is a process which is not only more straightforward but 
also considerably more economical than calculating the energy as a function of distance and 
numerically differentiating it. A detailed description of the interaction between two identical 
particles using this approach will be presented elsewhere. 

* * * 
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